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Abstract 

An Inverse Scattering Method is developed for the Camassa-Hohn 
equation. As an illustration of our approach the solutions correspond- 
ing to the reflectionless potentials are constructed in terms of the scat- 
tering data. The main difference with respect to the standard Inverse 
Scattering Transform lies in the fact that we have a weighted spectral 
problem. We therefore have to develop different asymptotic expan- 
sions. 
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1 Introduction 

In this introductory section some well known facts about the Camassa-Holm 
(CH) equation and the related spectral problem will be highlighted. The 
CH equation 

Ut — Uxxt + 2a;iia; -|- 'iuUx — 2UxUxx — UUxxx = 0, (1) 

where a; is a real constant, gained popularity as a model describing the 
unidirectional propagation of shallow water waves over a flat bottom [511251 

^On leave from the Institute for Nuclear Research and Nuclear Energy, Bulgarian 
Academy of Sciences, Sofia, Bulgaria. 



1 



I27j as well as that of axially symmetric waves in a hyperelastic rod |16[ ll8j. 
It firstly appeared in [^Hl as an equation with a bi-hamiltonian structure. 
CH is a completely integrable equation E El [211 , describing 
permanent and breaking waves jlUl |S| . Its solitary waves are stable solitons 
if a; > jSlElEllSEl or peakons if a; = 0121 • CH arises also as an equation 
of the geodesic flow for the right-invariant metrics on the Bott-Virasoro 
group (if to > 0) and on the diffeomorphism group (if to = 0) [S^ ^1 |H1 ^1 . 
The bi-Hamiltonian form of (pQ) is [SI i2llj : 

mt = —{o — o) — — — - = —{2ujd + mo + am) — — -. (2) 
om dm 

where 

m = u — Uxx (3) 

and the Hamiltonians are 

i/i[rn-] = —J mudx (4) 

H2[m] = ^ j{u^ + uul + 2Lvu'^)dx. (5) 

The integration is from — oo to oo in the case of Schwartz class functions, 
and over one period in the periodic case. 

There exists an infinite sequence of conservation laws (multi-Hamiltonian 
structure) Hn[m], n = 0, ±1, ±2, . . ., such that [T9 l liH l IsHI [2^ 

-(d-o) — — - = -(2ujd + mo + dm) ; — - — -. (6) 

dm dm 

The equation (0) admits a Lax pair (HI 



^xx = + A(m + a;))^' (7) 

= (^-n)^, + ^^' + 7* (8) 

where 7 is an arbitrary constant. We will use this freedom for a proper 
normalization of the eigenfunctions. 

In our further considerations m will be a Schwartz class function, uj > 
and m{x, 0) + lo > 0. Then m{x, i) + w > for alH |0]. For a discussion of 
the periodic case we refer to ^3] and |7|. Let k"^ = — Xu), i.e. 
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The spectrum of the problem under these conditions is described in 
[5]. The continuous spectrum in terms of k corresponds to /c - real. The 
discrete spectrum (in the upper half plane) consists of finitely many points 
kn = inn, n = 1, . . . , N where k„ is real and < k„ < 1/2. 

For all real k ^ a basis in the space of solutions of can be introduced, 
fixed by its asymptotic when x — > oo 

iPi{x,k) = e-*^^ + o(l), x^oo; (10) 
Mx,k) = e'''^ + o{l), x^oo. (11) 

Another basis can be introduced, fixed by its asymptotic when x — oo: 



ipi{x,k) = e-^'=^ + o(l), x^-oo; (12) 
ip2{x,k) = e*'=^' + o(l), x^-oo. (13) 

For all real /c 7^ if ^{x,k) is a solution of ((Tj), then ^{x,—k) is also a 
solution, thus 

(pi{x,k) = (f2{x,-k), il)i{x,k) = '4)2{x,-k). (14) 
Due to the reality of m in for any k we have 

Lpi{x,k) = i^2{x,k), tpi{x,k) = 4)2{x,k) (15) 

The vectors of each of the bases are a linear combination of the vectors 
of the other basis: 

ipi{x,k) = Tii{k)Tpi{x,k) (16) 
1=1,2 

where the matrix T{k) defined above is called the scattering matrix. For 
real k ^ 0, instead of ipi{x,k), ip2{x,k), 'ipi{x,k), ^2(2;, A;) due to (|T^ . for 
simplicity we can write correspondingly ip{x,k), ip{x,k), ilj{x,k), ^{x,k). 
Thus T{k) has the form 

- ( 'It! ) (^^) 

and clearly 

(fix, k) = a{k)'ip{x, k) + h{k)i}{x, k). (18) 
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The Wronskian W{fi, /2) = fidxf2 — f2dxfi of any pair of solutions of 
does not depend on x. Therefore 

W{ip{x, k), (p{x, k)) = Wiipix, k),ij{x, k)) = 2ik (19) 
From (UHl) and ^ it follows that 

|a(A:)p-|6(A:)|2 = l, (20) 

i.e. det(r(fc)) = 1. 

In analogy with the spectral problem for the KdV equation [^S], one 
can see that the quantities T{k) = a~^{k) and TZ{k) = b{k)/a{k) represent 
themselves the transmission and reflection coefficients respectively [HI 111). 
From ()2U() it follows that the scattering is unitary, i.e. 

|T(fc)|2 + |7^(A;)i2 = l. (21) 

The entire information about T{k) (|17() is provided by TZ{k) for A; > only 
|llj . It is sufficient to know TZ{k) only on the half line k > 0, since from 
JH and (HHI), a{k) = a{-k), h{k) = h{-k) and thus n{-k) = 1l{k). 

At the points of the discrete spectrum, a{k) has simple zeroes 0, there- 
fore and tp are linearly dependent (|18)) : 

99(x, IKn) = bnlpix, -lUn)- (22) 

In other words, the discrete spectrum is simple, there is only one (real) 
eigenfunction if'^"'\x), corresponding to each eigenvalue and we can 
take this eigenfunction to be 

^^''\x) = ^{x,iKn) (23) 

The asymptotic of according to (HU), ((221) is 

y,(")(x) = e''"^ + o(e'^"^), x^-oo- (24) 
(^W(x) = 6„e-'^"^ + o(e-^"^), x ^ oo. (25) 

The sign of 6„ obviously depends on the number of the zeroes of99("). Sup- 
pose that 

Q < Ki < K2 < ■ ■ ■ < KN < 1/2. (26) 
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Then from the oscillation theorem for the Sturm-Liouville problem H , V'^"'^ 
has exactly n — 1 zeroes. Therefore 

= {-ir-'\bn\. (27) 

The set 

S = {n{k) {k>0), Kn, \bn\, n = l,...N} (28) 

is called scattering data. The Hamiltonians for the CH equation in terms of 
the scattering data are presented in pHf] . 

The time evolution of the scattering data can be easily obtained as fol- 
lows. From (|18)) with x ^ oo one has 

ip{x, k) = a(A:)e-^^^ + b{k)e'^'' + o(l). (29) 
The substitution of ip{x, k) into (jSJ with a; ^ oo gives 



1 , , 

ft = + If (30) 

From (|29p . (|3flj) with the choice 7 = ik/2\ for the eigenfunction ip{x,k) we 
obtain 



a{k,t) = 0, (31) 

b(k,t) = jbik,t), (32) 

where the dot stands for derivative with respect to t. Thus 

a{k, t) = a{k, 0), b{k, t) = b{k, 0)e^*; (33) 



T{k, t) = T{k, 0), n{k, t) = n{k, 0)ex*. (34) 

In other words, a{k) is independent on t and will serve as a generating 
function of the conservation laws. 

The time evolution of the data on the discrete spectrum is found as 
follows, inn are zeroes of a{k), which does not depend on t, and therefore 
kn = 0. From (jHJ with 7 = ik/2X; k = iKn and (|25j) one can obtain 

-On- (35) 



1-4^2 
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The Poisson brackets for the scattering data of the Camassa-Holm equa- 
tion are computed in where also the action-angle variables are expressed 
in terms of the scattering data. 

In Section 2 we compute the asymptotics for large k of the scattering 
data and the eigenfunctions, which we use in Section 3 to develop the In- 
verse Scattering Transform for the CH equation. A number of recent papers 
[HI Uni 1^ OH used a Liouville transformation to reduce the weighted spec- 
tral problem Q to a standard problem. Our approach is more direct, and 
provides, we believe, more transparent formulas. The special case of reflec- 
tionless potentials {T^{k) = for all k) which corresponds to the important 
class of solutions, namely the multi-soliton solutions is separately studied in 
Section 4. A formula for the A^-soliton solution is obtained. 



2 Analytic solutions and Riemann-Hilbert Prob- 
lem 

For the application of the Inverse Scattering Method it will be necessary the 
asymptotics for large k of a{k) and the Jost solutions to be found. Firstly 
we compute the asymptotic of a{k). 

The solution of ((T)) can be represented in the form 

ip{x,k) =exp(-ikx+ xiy,k)dy]. (36) 

J —oo 

For Im A; > and x oo, ip{x, k)e^^^ = a{k), i.e. 

POO 

lna(A;) = / x{x,k)dx, Im A; > 0. (37) 



Since a{k) does not depend on t, the expressions J^^x{x,k)dx represent 
integrals of motion for all k. The equation for x{x,k) follows from (jT)) and 

dsni) 

xAx, k)+x^- 2ikx = (k^ + \)m{x) (38) 
and admits a solution with the asymptotic expansion 

oo 

P-n 



X{x,k) =pik+po + ^^. (39) 



n=l 



The substitution of (IHU]) into gives the following quadratic equation for 
Pi- 

pI - 2ipi + — = 0, (40) 

U) 
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with solutions 



pi=i(l±,/l + ^) 



(41) 



Since f^^Pi{x)dx is an integral of the CH equation, presumably finite, we 
take the minus sign in (|41() . One can easily see that po and all p-2n are total 
derivatives |24j and thus we have the expansion 



lna(A;) = —iak + 



n=l 



(42) 



where a is a positive constant (integral of motion): 



a 



m(x) . , 
1 + — ^ -l]dx, 



UJ 



(43) 



and /_„ = f^^P-n are the other integrals, whose densities, p_„ can be 
obtained reccurently from ^ [^ITT]. 

In terms of the scattering data a can be expressed as 



TV 



^ \1-2k. 

n=l 



ln|a(fc)| 
71" 7o 4A;2 + 1 



dk. 



(44) 



The asymptotic of a{k) for Im /c > and |A;| oo from (|42() is a(A;) 



e or 



e*"*=a(A:) ^1, Im A: > 0, |A:| ^ cx). (45) 
When A; is in the upper half plane the following expression is valid 



N 



1 



In a{k) = —iak + In r— H : 



n=l 



^"""(^'HdA:'. (46) 



Let us now consider the asymptotic of the Jost solutions, starting for 
example from ip{x, k), ()1U() . One can check that the asymtotic for oo 
has the form 



tp{x, k) 



-ikx+kG(x) 



r]{x, k) 



/ , N / N -^^1(2;) X2(X) 

ri{x,k) = Xo{x) + ^^ + ^)^ + ..., 



k 



k^ 



(47) 



where, due to (fT(l]) . G{x) and r]{x, k) ^ 1 for x ^ 00. The substitution 
of (EZI) into gives explicitly G{x), Xq, Xi, . . .: 



■ip{x 



' L\?7l(x) + Uj) 



UJ \l/4 Xiix) 

+ ^ + 



(48) 
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Introducing the function 



^(x) = exp 



, m(y) +UJ , ■ 
x+ I l)dy 



to 



(49) 



which looks hke a deformation of the ordinary exponent, H48|) can be written 
as 



-ik 



e(x)v/2^^iM + ^M + ... 



i'{x) 



k 



A;2 



(50) 



Furthermore, the function x{x, k) = '([^{x, k)e^^^ is analytic for Im A; < 0, 
[S]. This follows from the representation 



x{x, k) = 1 — — / : 'm{x')x{x' , k)dx' . 



(51) 



Notice that 



X ( m.{y)+u) 



1 is bounded for all values of x. Indeed, 



m{y) + uj 



Udy 



m{y)dy 



^ uj[l + 



< 



\m{y)\ 



UJ 



since 'm(x) is a Schwartz class function. Therefore the function 

i^{x,k) = i;{x,k)[^{x)y'' 

is also analytic for Im k < 0. 
Similarly, 



dy < oo 



(52) 



ip{x,k) = (/?(x, /c) exp I 



ik 



X + 



■m{y) + bj 



UJ 



l]dy 



e(x)v/2^^iM + ^M + ... 



e(x) 



k 



k^ 



(53) 



is analytic for Im A; > 0. 

Multiplying ^ by (,{x)/a{k) and using (gHl), (ES) and ((SSI) ^e obtain 



(p{x,k) 



e*^°a(fc) - 



ip{x,k) +n{k)ip{x,k)[£,{x)] 



2ik 



(54) 



The function ip{x, A;)/(e*^"a(A;)) is analytic for Im A; > 0, ip{x, k) is analytic 
for Im < 0. Thus, (|54j) represents an additive Riemann-Hilbert Problem 
with a jump on the real line, given by TZ{k)'i/j{x, k)[^[x)]'^^^ . 
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3 Integration of the CH equation by the Inverse 
Scattering Method 

Let us take an arbitrary k from the lower half plane (Im A; < 0). Then using 
the Residue Theorem, and (|^^ we can compute the integral 

1 f ^{x,k') dk' _ ip{x,iKn) 
2vri Jc^ e^^''^a{k') k' - k ~ ^ - A;)e-'^""a'(iK„) 

^ a'{iKn){iKn - k) 



n=l 



where C+ is the closed contour in the upper half plane (Fig.l). On the other 
hand, due to (|5H) the same integral can be computed directly as 
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1 / ip{x,k') dk' 



2vri Jc, e*'' "a(/c') k' - k 



1 1"^ / _ ,\ dk' 

' A;') + n{k')ijix, k') [^ix)f 



1 r ^{x, k') dk' 

2m e^'^'^aik') k' - k' ^ ' 

where r+ is the infinite semicircle in the upper half plane (Fig.l). Using the 
expansion (|53() and (|45l) . it is straightforward to compute that the integral 
over r+ is simply (l/2)(^(a;)/^'(x))^/^ 
Similarly, 



k 

= 11 + ^ /_ 

where C_ is the closed contour in the lower half plane, r_ is the infinite 
semicircle in the lower half plane (Fig.l). Due to (|5U() . ()52() the integral over 

r_ is (i/2)(e(x)/e'(x))V2. 

Now, from ((HSJ) ~ it follows that for Im k < 0, 



^ a'{iKn){k - iKn) 
n=l 

The expression ((55|) . taken at /c = — ZKp, p = 1, . . . , gives 



= (1^) + / n{k')l{xXmx)] 



^^'[x)J J-oo k' -\- IKp 

.^ 6n[e(x)]-2'-"V>(x,-m„) 
+*2^ 77^ — \( ~_L — ^ ■ 

The equations 1)581) - (|59() represent a linear system, from which 'ijj{x,k) 
(for real A:) and V'(^) —iHn) can be expressed through which, indeed is yet 
an unknown function. 

Let us now recall that \{—i/2) = 0. Since iIj{x, k) does not depend on m 
for A = and since ^'(a^j k) is defined by its asymptotics at — oo, it follows 
that ^(x, -i/2) = Thus, for k = -i/2, ^ gives 
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a'(iAc„)(K„ + l/2) ■ ^ ^ 

Since ip{x,k) and —ii^n) are known from (|58|) - (|59|) . the equation 
()6Up is a first order differential equation for ^, wliich can be directly inte- 
grated to give £,{x). In other words, ([58]) - (|6fl|) represent a system of singular 
integral equations for ijj(x,k), ij{x, —iHn) and ^(x). 

Since the time evolution of the scattering data is known (|.S5)) . the depen- 
dence on i.e. ^(x, t) is also known, expressed by the scattering data. Thus 
the set S, p8|) uniquely determines the potential: from (j35|l one obtains 



m(x, t) = ijj 



e(x,t) 

4 Reflect ionless potentials 



(61) 



The inverse scattering is simplified in the important case of the so-called 
reflectionless potentials, when the scattering data is confined to the case 
TZ{k) = for all real k. This class of potentials corresponds to the A^- 
soliton solutions of the CH equation. In this case b{k) = and \a{k)\ = 1 
H21|) and from ()46|1. (|44() one can easily find that ia'{inp) is real: 



where 



a 



2^'"' n 



E."({ 



1 + 2k„\2 



l^p -\- Kyi 



n=l 



Thus, ia'{iKp) has the same sign as 6„, (|27)) and therefore 



> 0. 



ia'{iKp) 

The time evolution of c„ due to is 



(62) 



(63) 



(64) 



Cn{t) = Cn (O)exp 



1-4^2 



(65) 



The equation (|59|) represents a linear system of equations for the quan- 
tities V'(x, —iKp): 
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N 



tpix, -iKp) + ^ 



■^{x,-iKn) = (jri-T) > p= l,...,iV. (66) 



n=l 



or 



g(x) N^l/2p 



where 



Kp -\- Kyi 



i.e., finally 




(67) 



(68) 



V'(x, -iKn,t) 



N 



Now (jnni) gives (note that from ^ ^{-oo,t) = 0) 



Y,A~p[C,tl n = l,...,N. (69) 
p=i 



x = x(e,t) = in (i-Yl 



cn(t)e 



(70) 



the time evolution of c„ is known (|65j) . This represents an implicit relation 
from which can be expressed as a function of x and t. Thus, in this case the 
scattering data uniquely determine ^ = ^(x, t) and therefore the potential 

m{x,t) (jnu. 

Using (|70j) and (|61j) we obtain the A^-soliton solution. Indeed, for fixed 
coordinates xq, to, since x is a monotonically increasing function of ^ there 
is a unique > (which is treated as a parameter from now on), such that 
Xq = X(^O)io)- Furthermore, we have 



(71) 



and from here and (|6H1 



m{xo,to)=m{X{^o,tQ),to)=uj [CoX^{Co,to)) -1 



(72) 
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u{xo,to) = u{X{^o,to),to) 
Finally, the A^-soliton solution is 



uix,t) = j\~\-~^^^^^)\c^X^\i,t)dS,-u:. (73) 

Note that is an explicitly defined function ((70)) in terms of the scat- 

tering data. Thus, the solution (f73|) does not depend on any additional 
parameter. 

For example, for the one-soliton solution the function X{^,t) is: 



=ln 



■2ki 



1 

Kl + 1/2 



-2ki 



(74) 



Note that since ^ - ^r^q^ > [cf. ^] and ci(t) > [cf. (jni], both the 
nominator and the denominator in ()74j) are always positive and singularities 
do not appear. 



5 Conclusions 

In this paper the Inverse Scattering Method for the CH equation is outlined 
in the case when the solutions are confined to be functions in the Schwartz 
class, a; > and 771(2;, 0) + a; > 0. The A^-soliton solution is explicitly 
constructed. The inverse scattering based on a Liouville transform, which 
leads to a spectral problem of the standard form —^yy + Q{y)^ = fj.^ (the 
Schrodinger equation) is developed in series of works [HI IMl IS ES] • In this 
case Q and m are related through the Ermakov-Pinney ordinary differential 
equation. The construction of the soliton solutions based on the bilinear rep- 
resentation of the CH equation (Hirota's method) is presented in 36, 37.. 55], 
see also [HS;, for a parametric representation of the A^-soliton solution. The 
situation when the condition m{x, 0) + a; > on the initial data does not 
hold is more complicated and requires separate analysis j2HlE|- lirn{x, 0)+uj 
changes sign there are infinitely many positive eigenvalues accumulating at 
infinity, and singularities can appear in finite time in the form of wave break- 
ing (inf^gnjita;} — > —00, while u stays uniformly bounded), cf. [01 Ell- The 
inverse scattering for multi-peakon solutions (for their existence we must 
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have w = 0) is reported in ^[2]. For the periodic solutions see (in 
this setting a scaling transform shows that there are no qualitative differ- 
ences between the cases uj = and w 7^ 0). The traveling- wave solutions of 
the CH equation (w = 0) are classified in j39j . The Darboux transform for 
the CH equation is obtained in 41 . The construction of multi-soliton and 
multi-positon solutions using the Darboux/Backlund transform is presented 
in [111123]. 
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